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Abstract 

We propose a class of forward-backward asymmetries with respect to the subprocess 
cm. scattering angle in V+jet production at hadron colliders, with V being any of 
(W, Z, 7, H), which are directly proportional to the gluon distribution g(x). The infor- 
mations that these asymmetries can provide are complementary of those reachable from 
measurements of the transverse momentum and/or the rapidity distributions of a V in 
kinematical regimes where gluon scattering subprocesses dominate. The accuracy which 
can be reached in the W, Z and 7 cases, at the upgraded Tevatron and at LHC, should 
allow a considerable improvement of our knowledge of the gluon distribution function, 
especially at large x. 
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1 Introduction 



One difficulty in our effort to test the Standard Model (SM) at a hadronic collider, and 
search for any new physics (NP) beyond it, is due to the uncertainties pertaining to 
the gluon distribution |l|, |2|, [3|. Remembering for example, that gg —>■ H provides the 
dominant contribution to H-production at LHC in the region of 100 GeV < nin < 
800 GeV 0, |J, we infer that a good knowledge of the gluon distribution is necessary for 
studying the properties of the Higgs particle and precisely estimating the backgrounds. 

The main constraints to the gluon distribution inside a nucleon at present, arise from 
DIS measurements which probe the low x range ||, and at higher x from measurements 
of the pt distribution of a prompt photon produced in pp — > 7X at ^/s = 23 GeV WA70 
0, and in pBe — > 7X at pi a b = 530 GeV |§, and also from heavy flavor production [|]. 
The combinations of subprocess cross sections contributing to prompt photon production 
are always forward-backward symmetrical in the subprocess cm., and consist of 

ddjqq -» jg) m 

di ' 1 1 

and 

do{gq{q) ~> jq(q)) + da(q(q)g -> jq(q)) ,^ 
dt dt 
as well as of subprocess cross sections in which the photon comes from a quark or gluon 
fragmentation. The sensitivity of such measurements on g(x), arises from the fact that 
for pp or proton-Nucleus scattering, the second subprocess (0) dominates over most of the 
Pt region [[J , provided that the produced photon is constrained to be sufficiently isolated 
JTTJ, in order to reduce the magnitude of the fragmentation contribution. 

Further information on g(x), from such forward-backward symmetrical subprocesses, 
could arise at LHC and the upgraded Tevatron. It has been shown [1 l.| that a detailed 



measurement of the rapidity distributions of an inclusively produced gauge boson W, Z or 
7, accompanied by a hard jet might provide an accurate determination of the gluon struc- 
ture function. In principle, the same kind of information could also arise from Higgs+jet 
production, but the expected statistics is very limited in this case. Since the theoretical 
treatment of any (V+jet) production for (V = W, Z, 7, H) is very similar, we consider 
below all these cases together. 

Thus, if in a (V+jet) pair production, we restrict ourselves to measurements of the 
Pt and/or rapidity distributions of V, then we are only sensitive to the combinations of 
subprocess cross sections^] 

da(qq -> V g) 

di u 

and 

dv(gq(q) -> Vq(q)) + da(q(q)g -> Vq(q)) ^ 
dt dt 



1 For the moment we disregard contribution from photon fragmentation, to which we come back below. 
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while for Higgs production we also get contributions from 



ddjgg -* Hg) 

di u 

These combinations of the subprocess cross sections lead to (V+jet) distributions which 
are symmetrical with respect to the (V+jet) cm. rapidity y, as well as with respect to 
cosd*, where 8* is the subprocess cm. scattering angle. 

We want to stress in this paper that in addition to these quantities, there exist a 
contribution proportional to the combination 

dv(gq(q) ~> Vq(q)) _ da(q(q)g -> Vq{q)) 
dt dt 

of the subprocess cross sections, which induces (V+jet) distributions which are anti- 
symmetric with respect to y, as well as with respect to cos#*. It turns out that these 
distributions are directly sensitive to g(x) (in fact they are directly proportional to it) 
and supply independent additional information, especially in the large x range, which to 
our knowledge has not been used so far, certainly due to the lack of statistics in present 
colliders. Of course, they also provide new information for the other parton distributions. 
The aim of the present paper is to study these quantities. To achieve this goal we con- 
struct for each (V+jet) final state, a forward-backward asymmetry with respect to cos#*, 
which is thus a function of the subprocess energy squared s = M 2 , as well as an odd 
function of the rapidity y of the center of mass. 

In calculating the aforementioned asymmetries for the V = W, Z, H cases, our 
philosophy is to keep the leading QCD contribution for the (V+jet) final state, and 
include the antenna pattern effect arising from the additional (V+jet+ soft gluon) final 



state integrated in a suitable phase space region |12], |13|], as a rough estimate of higher 
order QCD corrections. It turns out that the exact size of this antenna phase space 
region, is not very important for the considered asymmetries. For the (7+jet) case, 
we also include the photon fragmentation contribution, whose effect is again not very 
important, provided that a sufficiently strong isolation cut is imposed on the produced 
photon. Such an isolation is anyway desired, in order to increase the sensitivity of the 
above asymmetry on the gluon distribution. 

The content of the paper is the following. In Section 2, the formalism is presented 
containing the definition of the forward-backward asymmetry for the V+jet production, 
and including also discussions of the possible antenna pattern effects and the photon 
fragmentation contribution. In Sect. 3 the sensitivity of the above asymmetry to the 
gluon distribution is discussed, while the conclusions are given in Sect. 4. Finally the 
subprocess cross sections and quark distributions for the various V cases, are given in the 
Appendix. 



2 Formalism 
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2.1 The V+jet contribution. 

The generic subprocess contributing to the pp — > V jet ... (V = W, Z, 7, H) cross section 
(apart from the photon fragmentation case discussed separately), is written as 

a(pi) + 6(p 2 ) -> V(p 3 ) + c(p 4 ) , (7) 

where the momenta are indicated in parentheses and the masses of the partons (a, b, c) 
are neglected. As usual s — {p\ + P2) 2 , t = (p^ — Pi) 2 , u = (p^ — p?) 2 . We also define 
r = s/s, where s is the LHC (or Tevatron) cm. energy- squared. The rapidity y of 
the cm. of the (V+jet) subprocess determines the momentum fractions of the incoming 
partons through 

x a = y/re y , x b = y/re~ y , (8) 
while the scattering angle in the cm. of the same subprocess satisfies^] 

cos r = in* . (9) 

u + t 

The (V+jet) production cross section in pp collisions is given by 
daipp - V jet) s - ml f Q2) da{gg _ V,) 



drdyd cos #* 2 [ rft 

+ Sy(x a ,X 6 )— CT ^ / ' / 



eft 

+2 fo^a, Q 2 )Si/(x 6 ) + Y, v (x a )g{x b , Q 2 )] 

+^ fo^a, Q 2 )£y(x 6 ) - Sy(x a )5f(x 6 ,<5 2 )] 



cicr(^g -> Vq) da(qg Vg) 



eft 

eft>(#<? -> Vg) do-(e/p -> Vg) 



(10) 



where ([8], |^) should be used. In flip]) , Q — Pr/2 is the usually preferred factorization 
scale of the distribution functions which mimic the next to leading order corrections. The 
Sy,Sy terms describe combinations of quark (antiquark) distributions weighted by the 
quark electromagnetic or weak charge, while da denote the correspondingly normalized 
subprocess cross sections, for the various V — H, W, Z, 7 cases. These are given in 
the Appendix for a pp Collider, while for the pp Tevatron they should be modified in an 
obvious way. 

The last two terms of eq.(|l(]) have been constructed by combining the corresponding 
gluonx quark and the quarkxgluon contributions in a symmetrical and an antisymmetrical 
parts. It is important to note that the last term in ( ]T0| ) is antisymmetric with respect 
to y or cos 8* for the Tevatron or LHC Colliders, while all other terms are symmetric 
in both these variables. This last term will therefore be washed out if we integrate over 
y, as it always happens whenever we only look at the rapidity and/or px distribution 

2 For a collection of relevant kinematical formulae see e.g. Appendix A3 in [p^| . 
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of a single V or jet. It will also be washed out, in case we cannot discriminate V from 
the accompanying jet, and therefore also in dijet production. Particularly for the dijet 
case, we note that so long we cannot discriminate between a gluon and a quark jet, the 
cos ^-antisymmetric term in (|10D vanishes identically, and only the symmetric part of 
([Top contributes. For studying therefore the last term in (|T0|), which is the main purpose 
of the present paper, the two high p T objects in the final state should be distinguishable, 
as e.g. in the case &W, Z or 7 production, accompanied by a high p T jet. 

We next turn to the corrections to the various terms in (|T0|), arising from soft gluon 
emission in the antenna approximation, and from photon fragmentation in case V = 7. 



2.2 Adding the antenna pattern contribution. 



The hadronic antenna patterns in V = W, Z, 7 + jet production have been shown 
to provide a valuable diagnostic tool for probing the nature of the underlying parton 
subprocess. They could also provide a good tool for distinguishing between conventional 
QCD and new physics production |15| (large Et jet events in hadronic collisions, or 
large anomalous Q 2 events at HERA). Here the philosophy is different: in order to make 
our treatment more realistic, we add to the leading order calculation of the V+jet cross 
section, the contribution from the soft gluon emission. Thus, for each subprocess like in 
(|7|), we consider the corresponding subprocess 



a(pi) + b(p 2 ) -> V(p 3 ) + c(p 4 ) + g(k) 



The soft gluon emission in such subprocesses is controlled by the basic antenna pattern 
distribution [|12| 

n = p i Pj , , (12) 

Pi ■ k pj ■ k 

which should be adequate for the description of gluonic minijets with energies fco much 
smaller than those of the hard jets [|l^ . 

The contribution from the (y+jet + soft gluon) production cross section in pp colli- 
sions is then given by (compare ([T0|)) 



da{pp — > V jet) 
drdyd cos 9*d? k 



3a s (s 



8tt 2 



^(x a ,Q 2 )^(x b ,Q 2 )([12] + [14] + [24]) 



da{gg -> Vg) 



dt 



+ E v (x a ,x b )[[U] 



[24 ] _ M] da ^' 



Vg) 
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dt 



1 



+ 2 [d( x a, Q 2 )^v(x b ) + Y, v (x a )g(x bl Q 2 )] 



[12] + [14] 



[24]\ d&jgq -> Vg) 
9 J dt 



F [12] + [24] 



[14]\ dajqg -> Vg) 

9 ) di 



+ 77 [9(? a , Q 2 )S y (x 6 ) - i: v (x a )g(x b , Q' 
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[12] + [14] _ m *><*> r m _( [12]+ m _ m r w 

9 / dt V 9 J dt 



.(13) 



Neglecting all parton masses, except my, we use the notation p± = V s/2(l, 0, 0, 1), 
= Vs/2(1, 0, 0, —1) for the description of the momenta of the incoming partons a, 6 in 
their cm. (compare flllD ), and the notation 



?2 



l-i 

Pv 

Pi 



Pz = (E VT coshy v , p T , 0, E VT smhyy) , 
pr(cosh^*, —1, 0, sinhr/*) , 

/c r (cosh(r/* + 5rj), cos(S(p), sm(5(f), sinh(?7* + 5r])) , 



(14) 
(15) 
(16) 



for the final state particles in the same frame [ 13(1 . Here Eyr = a/ my + gives the 
transverse energy of V, and (y v , rj*) denote the rapidities of (V, c) in the (a, 6)-c.m. 
frame. The final gluon is taken to be soft, which means 

(17) 



Aqcz? < h < vs/2, EyTCOshy^, p T cosh?7* 



Note from (|15|, (T6|) , that to the extent that 5r] , 5ip are small, the rapidity and azimuthal 
angle of the minijet generated by the soft gluon are close to those of the hard jet generated 
by parton c; compare (|ll]). 



In terms of these variables the antenna pattern coefficients in (|i~2| ) may be expressed 



as 



[12] 
[14] 
[24] 



,,5il 



k\ [cosh(<5?7) + cos(5v?)] 

k\ [cosh(5^) + cos(5<£>)] 
while the phase space of the soft gluon is determined by 

d 3 k 

— — = krdkr dSrj d5ip . 
k 



(18) 
(19) 
(20) 



(21) 



Defining 5R = a/ 5t] 2 + 5(p 2 , we choose to integrate the antenna activity contained in 
( |T3"D in the region 



< 5R < 5R 2 , 

kTmin — — ^Tmax ■ 

Adding then ([TOD and the result of integrating (|l"3"D we get 



(22) 
(23) 



da{pp — > V jet) s — m v 
drdyd cos 9* 2 



g(x a ,Q 2 )g(x b ,Q 2 )I, 



99 



dvj gg -> Vg) 
di 



E v (x a ,x b )I, 



da(qq' -» Vg) 
dt 



+ 2 Q 2 )^v(xb) + E v (x a )g(x b , Q 2 )] I, 

+ 2 W 2 ^' Q 2 ) e v(e&) - S v (x a )5((x fe , Q 2 )] I, 



m 



gg 



da(gq -> Kg| + df(Q9^ Vq) 



dt dt 
da(gq -> Kg) ^(gg -> Kg) 



where the antenna contribution is contained in the parameters 

5 \ 



' gg 



L gg 



3a s ( k^ max 

1 + ^^rin 



1 + 
1 + 



2tt 2 
30^ 

3a 



In 



■In 



kTmax 
kTmin 
^Tmax 



27T 2 "* V k Tmin 
expressed in terms of the function| 



[ix(5R 2 2 -5Rl)+£(5R 2 , 5R } 



ir{5R 2 2 - 5R\) + -£(5R 2 , SRx) 
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{8R\ - 8R{) + £(5R 2 , 5Rx) 



Z(5R 2 , 6Ri 



d5R 5R 



8R-, 



2tt 



cosh(5i? cos( 



cosh (£i? cos 0) + cos(5i?sin. 



(24) 



(25) 
(26) 
(27) 



A general comment is in order here. The quantitative predictions on the magnitude of 
the colour flow are based on the hypothesis of Local Parton Hadron Duality jL6|, which 
assumes that colour coherence effects survive the hadronization stage; an hypothesis well 



confirmed by existing data [17 



2.3 Photon fragmentation contribution. 



The result ([24]) should be adequate for the V = W, Z, H cases. For the 7+jet case though, 
we should add the contribution from processes where the final state photon comes from 
the fragmentation of a quark or gluon. Although such processes are formally of higher 
order in QCD, they may occasionally turn out to be quite important, since the logarithmic 
growth of the fragmentation function due to scaling violations, compensates one power of 
a s |T2|]. The contribution from such fragmentation to the 7+jet cross section arises from 



each partonic subprocess 

a(Pi) + b(p 2 ) -> c(p 3 ) + d(p 4 ) , (29) 

in which the parton c fragments subsequently to a photon, through a fragmentation 
described by the function DJ(z,Q 2 ), where z denotes the fraction of the c-momentum 

3 Notice that this result depends crucially on the fact that the antenna integration region (f2|, ^3|) is 
chosen so that 5<p in (|l^) remains small, so that fc M never goes so close to p%, to make the perturbative 
treatment unreliable; compare (flq, na). 
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carried by the photon. Following [10] we write the next to leading QCD order calculation 
of the fragmentation contribution to pp — > 7 jet as 



daf ra9 {pp — * 7 jet) 



drdyd cos 8* 



^2 fa/p( X *' Q 2 )fb/p(%b, Q 2 



abed 



da(ab — > cd) 
di 



dzD2(z,Q 2 f ) 



(30) 



where f a / p (x a , Q 2 ) denotes the a-parton distribution function in a proton, at a scale Q ~ 
Pt/2. Yliabcd re f ers to the summation over the full list of subprocesses, q(q) + q(q) — > 
q(q) + q(q), q(q) + g — > q(q) + g, qq — > — > The expressions for the subprocess 
cross sections can be found in |Tj| |T2"|| . 

In (p0[), Q/ denotes the scale of the fragmentation function which, according to the 



next to leading order calculation in ||10|| , is determined by 



Qf 



pIr 

cosh ?7 7 



(31) 



where 7]^ is the photon rapidity in the laboratory frame, and p\ its transverse momentum. 
In (0), R = ^{5ri) 2 + (Sip) 2 denotes the size of the isolation cone around the photon 
produced from the c-fragmentation, within which the hadronic energy Ehad is constrained 



to be smaller than some max(Ef la d), which in turn determines also z m i n in (p0[ ) through 



max(E had ) 
E, 



(32) 



The smaller m.ax.(Eh a( i) is chosen, the more isolated the photon becomes, which in 
turn means a smaller photon fragmentation contribution to the cross sections. In the 
calculations below we use R = 0.7 and max(Ehad) — 4 GeV f2(| in (31 , [30|), which means 



constraining the photon to be quite isolated. Using then standard leading order expres- 
sions for the various partonic cross sections in (|30|) fl2]| , and the photon fragmentation fit 



zDJ{z,Q)) 



fx 

zD](z,Q 2 f ) = - 



2.21 - 1.28^ + 1.29^ 



v 0.049 



+0.002(1 



1.631n(l-z) 



z) 2 z~ x - u 



a 
2^ 



0.0243(1 



1.03^-0.97 



HQ 2 /K 2 QCD ) 



(33) 
(34) 



(with Aq CD = O.OAGeV 2 ), we calculate the photon fragmentation contribution to the 
cross sections, which should be added to the cross section in (^). It turns out that for 
all numerical applications presented below, da^ ra9 (pp —>■ 7 jet) from (|30"D is always less 
than 10% of da(pp — > 7 jet) from (f2"4|). 
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We are aware that our estimate of the fragmentation contribution is very approxima- 
tive. Nevertheless there are two basic reasons to be convinced that we are not far from a 
precise next to leading order calculation. Firstly, concerning photon fragmentation func- 
tions, it has been shown in |E| that, in the large z domain we probe (z ~ 1), the photon 
fragmentation functions beyond leading order do not strongly deviate from the leading 
log parametrization we have used. Secondly, with respect to the choice of the scale of the 
fragmentation function and its relation to the isolation cone, we should admit that a fully 
satisfactory treatment is lacking at present [23|. In the present work we simply followed 
the results of iTO . 



2.4 The Asymmetry 

The forward-backward asymmetry is defined as 



A v (r,y) 



J 6 d cos 8* 


dcr(pp^V jet) 
drdyd cos 9* 


dcr(pp^V jet) 
drdyd cos 9* 

9* 


ir-9*- 


6 d cos 6* 


da(pp—*V jet) 


. da(pp^V jet) 


TT-9*. 


drdyd cos 9* 


drdyd cos 9* 

9* 



(35) 



where e is a small positive number serving to exclude the angular region around the beam 
direction, where the cross sections are not measurable, and the pertubative treatment 
not applicable. The V production cross section appearing in fl3"5| ) is given in ([TU|) , and 
can be refined by including in it the gluon bremsstrahlung contribution at the antenna 
approximation as indicated in (0), as well as the fragmentation contribution from ( j30[) 
for the photon case. 

As it can be seen from (24, 30), the asymmetry A v (r,y) is always an antisymmetric 
function of y. For the V = W, Z, H cases, where (fPj) contributes (containing the leading 
order and the soft gluon contributions inside the integrals /«), it can be expressed as 

A v (r,y) = 



[g(x a )Y, v (x b ) - T lV (x a )g(x b )} I gq J i 



v- 

9Q 



[g(x a )T, v (x b ) + E v (x a )g(x b )] I gq J^ q + + S y (x a , x b )I qq J^ q + g{x a )g{x b )I gg JV g 



(36) 



The antenna pattern parameters appearing here, are given in (p5| ~ pT|) , while those from 
the various subprocess cross sections are written as 



J, 



91 



J, 



m 



J, 



1 



l-e 



2 Jo 
1 ^ 



d cos 8* 



da(gq -> Vq) da(qg -> Vq) 



dt 



dt 



d cos 6* 



o 

l-e 



da(gq -> Vg) da(qg -> Vq) 



d cos 9* 



dt 

da(qq' -> Vg) 



dt 



dt 



(37) 
(38) 
(39) 



J. 
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d cos 9* 



(40) 



where of course ( |40|) is relevant only in the Higgs case. The y-antisymmetry of A v is then 
obvious on the basis of (|8|, |36]). Its direct dependence on the gluon distribution g(x) is 
also clear in fl36|) . 

It is possible to write a corresponding formula for the A 1 case, but this time the result 
is more complicated, since the fragmentation contribution involves many subprocesses. 



3 Sensitivity to the gluon distribution. 

In principle the above asymmetries supply additional independent information on the 
parton distributions, and should be used together with the other usual measurements 
to constrain these distributions. Since these asymmetries are proportional to the gluon 
distribution, we would like to explore here their usefulness for constraining g(x), assuming 
that the quark distributions have already been precisely determined by other means. Thus 
here, we are interested in using the above asymmetries for the various production processes 
in order to improve our knowledge on the gluon distribution g(x), particularly in the large 
x region, where most of the uncertainties lie. We proceed as follows. 

We first want to see what uncertainty our present ignorance of the parton distributions 
in general, and of the gluon distribution in particular, imply for the above asymmetry. 
This is done by using the existing fits for the parton distributions [1], @, [§, which provide 
an estimate of the bands inside which the asymmetry may lie, for various values of the 
V^+jet invariant mass M = and the rapidity y of the V+jet pair. These bands are 
given by the dotted lines in Figsjl], |], ||, § for the W ± , Z, 7 and H production cases 
respectively, while the full line determines the "averaged" asymmetry in each of these 
cases. In all figures we use e = 0.05 to cut the forward or backward V production, with 
respect to the cm. angle; but the results are not sensitive to the exact magnitude of this 
value; (compare fl3"5|)). 

In the figures we also include the antenna pattern contribution generated by a soft 
gluon jet "close" to the hard jet generated by the parton c; compare ([16|, |22|, |23[). This is 
integrated in the phase space region [13|, |24fl : 



0-7 < 6R < 1 , 70GeV < k T < 200GW for LHC (41) 
0-7 < 5R < 1.3 , WGeV < k T < 30GeV for Tevatron (y£ = 2TeV) . (42) 

The effect of such antenna contributions is to increase the forward and backward cross 
sections and the asymmetry by roughly 10%, in all cases. 

For the 7 production, we also include the fragmentation contribution from (|30|). It 
turns out that for the strong photon isolation constrain imposed by ma.x(Eh a d) = ^GeV 
(compare (^, |32|), the fragmentation contribution to the forward or backward cross sec- 
tions are at the 10% level. This photon fragmentation contribution tends to decrease the 
asymmetry by roughly 10%. 
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Thus, the bands in Figs. 1-4, give an estimate of the present uncertainty on these 
asymmetries, which is due to our present ignorance of the gluon distributions. To see what 
can be achieved by LHC and/or the upgraded Tevatron, we indicate on the same figures 
the uncertainties of a possible (future) measurement of these asymmetries, assumed to lie 
on the "averaged" solid line. These uncertainties include the statistical ones implied by 
the integrated luminosities (2 experiments with 3 years of running) of and 12/fe -1 

for the LHC (y/s = 14 TeV) and the upgraded Tevatron (y/s = 2 TeV) respectively. For 
W and Z, we only retain leptonic decay modes with branching fractions Br = 0.16 and 
0.09, respectively. One should notice that the Z + jet final state is experimentally very 
clean. 

For 7 and H we assume a detection efficiency of 0.8 and 0.01 respectively. However 
7 + jet final state suffers from a background from neutral pions, which will limit the 
accuracy of such a final state. The detection efficiency in the Higgs case is of course very 
much dependent on the value of the Higgs mass which controls the value of the branching 
ratios in 77, WW*, ZZ* etc, and the value of 0.01 is only chosen for orientation. 

In the illustrations given in Figs. 1-4 we have computed the statistical uncertainties 
taking bins of Sy = ±0.02 (indicated by the horizontal lines in the figures) and 5r = 
(±300 GeV ■ M)/s. Eq.(|8|) gives the relation between the kinematical domains (r = 
M 2 /s, y) and (x a , Xj,). Note that reducing the size of the bins would make the domain 
in (x a , Xb) more restrictive, but it would simultaneously increase the statistical error. 

To these statistical uncertainties we quadratically add the uncertainties implied by the 
present knowledge of the quark distribution functions, estimated by their spread among 
the various models. One should take into account that at LHC the lepton pseudorapidities 
from weak gauge bosons decays will provide the key to measure the q, q distribution 
functions within an accuracy of ~ 1% [PTTj] . The resulting total uncertainties are indicated 
by the vertical error bars at the central point of each bin. 

Figs. 1-4 show the kinematical domains where the measurement of the asymmetries 
can provide useful constraints on g(x). In order to get informations on g(x) at large x, 
one needs measurements at large M and y; compare (§). The number of events is however 
decreasing with M and y. For LHC we see that asymmetries for W and 7 production are 
useful (i.e. at the percent level) for M < 1000 GeV and \y\ < 2, which means x values 
up to x ~ 0.3. The results for Z production are given in Fig.^. The expected sensitivity 
to g(x) is weaker than in the W, 7 cases and it can be useful only up to x ~ 0.1. The 
Higgs case is shown in Fig.[|. One cannot expect an improvement on the determination 
of g(x) from it, but it can contribute to consistency checks, particularly with respect to 
excluding possible New Physics contributions. 

For the upgraded Tevatron the lower energy (y/s = 2TeV) and the lower luminosity 
forces us to restrict to M < 200GeV, for which we could reach useful information (at the 
percent level) only for x < 0.1 and only from 7 production (compare Fig||), the W case 
giving only consistency checks. 
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4 Conclusions 



In this paper we have shown that forward-backward asymmetries with respect to the 
subprocess cm. scattering angle in V^+jet production at hadron colliders, are directly 
proportional to the gluon distribution function g(x). We have then studied how measure- 
ments of these asymmetries can improve our knowledge of this gluon distribution. 

We have considered V^+jet production at the upgraded Tevatron and at the LHC, for 
the cases V = W ± , Z, 7 and H. In order to make our treatment more realistic, we have 
added to the leading V+jet cross section, the contribution from the soft gluon emission in 
the form of an antenna pattern distribution. In the case of the photon we have also added 
the photon fragmentation contribution arising from the various subprocesses involving 
two body scattering among quarks and gluons. 

We have made an estimate of the accuracy at which these asymmetries can be mea- 
sured at the upgraded Tevatron and at the LHC in several bins of V^+jet invariant mass 
and cm. rapidity. Each bin corresponds to a certain domain in the fraction of momentum 
x carried by the partons. We have given the kinematical domains where this accuracy 
would allow a substantial improvement in the knowledge of the gluon distribution g(x). 
More explicitly: 

• At the upgraded Tevatron (y/s = 2 TeV, with an integrated luminosity of 12 fb^ 1 ) 
and using mainly the 7+jet process, one can expect to reach an accuracy of the 
percent level up to about x ~ 0.1. 

• At LHC (y^ = 14 TeV, with an integrated luminosity of 600 mainly through 
the 7+jet and VT+jet processes (and at a weaker level the Z+jet one) one should 
test g(x) at less than one percent up to x ~ 0.3. 

• A comparison of the results in the various modes, including the if +jet one (which 
is by itself not accurate enough to severely constrain g(x)) should give consistency 
checks and constrain possible non standard contributions. In this respect, the radi- 
ation antenna pattern may also help. 

The present study should just be taken as "a stone to be used for the construction 
of the building" . Obviously the experimental information provided by these asymmetries 
should be inserted in a global analysis involving all other observables and all parton 
distributions. In this context it should be remembered that the information contained in 
these asymmetries arises from subpro cess-terms which are completely different from those 
determining the present prompt photon |7|, || and the dijet production cross sections at the 



Tevatron ||25|| . The simple-minded analysis presented here was just aimed at showing the 
peculiar role of these asymmetries for the purpose of determining the gluon distribution 
function, provided - as it will be the case in future hadronic colliders - high statistics will 
allow to take into account this forward backward asymmetry in a global analysis. 
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Appendix: Parton expressions for V production at a pp Col- 
lider. 



For H production the relevant quark distribution to be used in ([10 , [24]) are 



E H (x) = £(x) = ]T[g(x,Q 2 ) + g(x,Q 2 )] , (A.l) 

q 

E H {x a ,x b ) = ^2 b(x a , Q 2 )q(x b , Q 2 ) + q(x a , Q 2 )q(x b , Q 2 )] , (A.2) 
i 

where the sum is over all light quark flavours. The subprocess cross sections can be found 



in 



19| , PH . Note that the gluon-gluon term in (|i~0D contributes only for Higgs production, 
and it should be ignored in the (W, Z, 7) cases below. 

For W production (W = W + + W~) 

E w (x) = E(x) , (A.3) 
T, w (x a ,x b ) = d(x a ,Q 2 )u(x b ,Q 2 ) +u(x a ,Q 2 )d(x b ,Q 2 ) + d(x a ,Q 2 )u(x b ,Q 2 ) 

+ u(x a , Q 2 )d(x b , Q 2 ) + [u^c , d -> s] , (A.4) 

provided we neglect the CKM matrix elements V t d, V ts . The corresponding subprocess 
cross sections, to leading order in QCD, are 

(A.5) 
(A.6) 
(A.7) 



da(qq' — 


Wg) 


a s V2G F M^ 


8 


t 2 + u 2 + 2M^s 


di 




4s 2 


9 


iii 


da(gq -> 


Wq') 


a s V2G F M^ 


1 


s 2 + u 2 + 2M^i 


di 




4s 2 


3 


—su 


da(qg -> 


Wq') 


a s V2G F M^ 


1 


s 2 + i 2 + 2M^u 


di 




As 2 


3 


-si 



Correspondingly, for Z production 

E z (x) = Y,(9l q + 9l q )[q{x ) Q 2 ) + q{x 1 Q 2 )\ , (A.8) 
q 

Z z (x a ,x b ) = ^2(g 2 Lq + g 2 Rq )[q(x a ,Q 2 )q(x b ,Q 2 ) +q(x a ,Q 2 )q(x b ,Q 2 )] , (A.9) 
q 

where the Z-couplings to the L and R quarks g^ q = tq — e^s 2 ^ and g Rq = —eqS^, 
are absorbed to the definition of the H z , distributions. In this normalization, the 
corresponding subprocess cross sections are given by (|A.5| - |A.7|) , after replacing M w — > 
M z . 

Finally for 7 production we have 

S 7 (x) = ^e 2 [g(x,Q 2 )+g(a;,Q 2 )] , (A.10) 
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^(xa.Xb) = ^2e 2 q [q(x a , Q 2 )q(x b ,Q 2 ) + q(x a ,Q 2 )q(x b , Q 2 )] . (A. 11) 

Since the square of the quark charge e 2 has been inserted in the definition (|A.10| , |A.11| ), 
it is removed from the subprocess cross sections, which is thus written to leading QCD 
order as 



da(qq — » 


7fl0 


7ia s a 


8 


i 2 + u 2 


di 




s 2 


9 


til 


da(gq — 


7?) 


na s a 


1 


s 2 + u 2 


di 




P 


3 


—sit 


dcr(qg — 


iq) 


na s a 


1 


s 2 + i 2 



dt 



-st 



(A.12) 
(A.13) 
(A.14) 
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Figure 1: Asymmetry for W^+jet production with invariant mass M, as a function of 
the pair rapidity y, for LHC (a,b,c) and the upgraded Tevatron (d). The band between 
the two dotted lines is due to the present uncertainty on the gluon distribution function. 
The solid line is the center of the band. The error bar shows the accuracy at which the 
asymmetry can be measured for each indicated bin (dashed lines) taking into account 
statistical errors and uncertainties on quark distribution functions. 
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Figure 2: Asymmetry for Z+jet production with invariant mass M, as a function of the 
pair rapidity y, for LHC (a,b,c) and the upgraded Tevatron (d). Same caption as in Fig.l. 
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Figure 3: Asymmetry for 7+jet production with invariant mass M, as a function of the 
pair rapidity y, for LHC (a,b,c) and the upgraded Tevatron (d). Same caption as in Fig.l. 
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Figure 4: Asymmetry for if+jet production with invariant mass M, as a function of the 
pair rapidity y, for LHC (a,b). Same caption as in Fig.l. 
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